Theory of integer quantum Hall effect in graphene 



o 
o 

5h . 

< 
OO 



13 
G 
O 

o 



> 

On 

m 
o 

O 

-i— > 
c3 



i 

c 

o 
o 



X 



r/m 



Igor F. Her but 

Department of Physics, Simon Fraser University, Burnaby, British Columbia, Canada V5A 1S6 

The observed quantization of the Hall conductivity in graphene at high magnetic fields is explained 
as being due to the dynamically generated spatial modulation of either the electron spin or the 
density, as decided by the details of Coulomb interaction on the scale of lattice constant. It is 
predicted that at a large in-plane component of the magnetic field such ordering may be present 
only at the filling factor / = ±1 and absent otherwise. Other experimental consequences of the 
theory are outlined. 



Graphene is a two-dimensional semi-metal with gap- 
less Dirac-like excitations near two points in the Bril- 
louin zone 0,0. When placed in the uniform magnetic 
field of several Tesla it exhibits plateaus in the Hall con- 
ductivity at integer filling factors / = An + 2, 0], [J], 
just as implied by the Landau level (LL) spectrum of 
the single-particle Dirac equation [Ej], with each LL be- 
ing four-fold degenerate in spin and sublattice indices. 
Although the strength of the Coulomb interaction be- 
tween the conducting electrons in graphene is similar to 
the bandwidth, the semi-metallic ground state at zero 
magnetic field indicates that it is below the critical value 
needed for insulation 0]. 

At magnetic fields above ~ 10T, however, additional 
quantum Hall (QH) states at / = 0, ±1,±4 appear ||. 
Plateaus in the Hall conductivity at other integers, such 
as at / = ±3 or / = ±5 for example, are at the same 
time conspicuously absent, even at the highest magnetic 
field of 45T. The experiment [1] also suggests that the 
activation gaps at / = ±4 are likely to be due to the 
Zeeman splitting of the first LL. The quantization at / = 
±1, however, implies that the fourfold degeneracy of the 
zeroth LL has been completely lifted, which calls for the 
inclusion of the Coulomb interactions into consideration. 
The theory of the integer QH effect in graphene would 
thus have at least to provide answers to the following 
questions: a) why has the fourfold degeneracy of only the 
n = LL been completely lifted at the magnetic fields 
and samples under study, b) why do new incompressible 
states at / = 0, ±1, ±4 require higher magnetic fields to 
appear than those at / = 4n + 2, and finally, c) what is 
the nature of the interacting ground states at different 
filling factors? 




The first question may be immediately answered by 
postulating that the Dirac fermions have acquired an ef- 
fective gap in a form of a 'relativistic mass' due to the 
Coulomb interaction (9j. Such a gap reduces the degen- 
eracy of only the zeroth LL. Including the Zeeman split- 
ting then leads to the spectrum of the effective single- 
particle Hamiltonian precisely as required by the ob- 
served pattern of quantization of Hall conductivity. It is 
unclear, however, under which circumstances, and which 
one of the multitude of such gaps that could exist, as 
discussed below, is generated dynamically. This problem 
is addressed here within the extended Hubbard model 
on a honeycomb lattice, with both on-site and nearest- 



FIG. 1: The proposed phase diagram of graphene in the mag- 
netic field. D z = g z B is the Zeeman energy, r is the chemical 
potential, and mo = 2g x B±/n is the characteristic size of 
the 'relativistic' many-body gap m. g x is the larger of the 
couplings in the CDW and AF channels (see the text). 



neighbor repulsions. This is the simplest Hamiltonian 
that mimics the effect of Coulomb repulsion and which 
contains the possibilities of charge-density-wave (CDW) 
and antiferromagnetic (AF) orders. Previous work indi- 
cated , that these two are in direct competition at 
zero magnetic field and when the interaction is strong. 
Here we solve the model in the physical weak coupling 
regime and in an external magnetic field, allowing only 
for the simplest AF state with the staggered magnetiza- 
tion parallel to the magnetic field. The results are sum- 
marized as the phase diagrams in Figures 1 and 2. At 
half-filling and for a weak enough Zeeman energy the 
system could be either a CDW or an AF, depending on 
which coupling dominates (Fig. 2). For a larger Zee- 
man energy the ground state at / = becomes magnetic 
with full lattice symmetry. Nevertheless, even in the lat- 
ter case increasing the chemical potential produces an 
incompressible state at / = 1, with the activation gap 
becoming equal to the 'relativistic' gap (Fig. 1). In con- 
trast, at weak Zeeman coupling we find a direct transi- 
tion between / = and / = 2 states. At / > 2 the 
'relativistic' gap vanishes. Experiments that would test 
the presented against other theories [12], 0, [3, 
are discussed. 
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FIG. 2: The phase diagram at half-filling (r = 0). The trans- 
lationally symmetric magnet exists for g c , g a < ng z B/2B±. 



We define the extended Hubbard model as H = Ho + 
Hi, where 

H = -t ui(^)v a {X+bi)+H.c., (1) 

A,i,a=± 

Hi = U ^n+(X)n_(X) + ^ n a (A)n cr i (A + 6^). 

(2) 

The sites A denote one triangular sublattice of the hon- 
eycomb lattice, generated by linear combinations of the 
basis vectors a x = (-s/3, — l)(a/2) and a 2 = (0, a). The 
second sublattice is then at B = A + b, with the vector b 
being either b x = (1/V3, l)(a/2), b 2 = (1/V3, -l)(a/2), 
or 6 3 = (— a/\/3, 0). a » 2.5A is the lattice spacing, and 
t w 2.5eV, C/ w 5 - 12eV, and E7/V « 2 - 3 

The spectrum of i?o becomes linear in the vicin- 
ity of the two Dirac points at ±K, with K = 
(1, l/V3)(2n/aV3) Q, [3- Retaining only the Fourier 
components near ±K one can write, in the continuum 
notation, 

Ho = f dSj^^K^holiDi^^x), (3) 

and 

7 (27ra) J 
«t (K + q),u\{-K + q),vU~K + q)), 

i = 1,2, where it was convenient to rotate the reference 
frame so that qi = q ■ K/K and q 2 = (K x q) x K/K 2 , 
and set Ti = e/c = vp = 1, where = is the 

Fermi velocity. Here 17071 = er 2 ® cr^, 17072 = —-^2 ® cr y , 
with ^2 as the 2x2 unit matrix, and a as the Pauli 



matrices. A w 1/a is the ultraviolet cutoff over which 
the linear approximation for the dispersion holds. The 
orbital effect of the magnetic field is included by defining 
Di = —idi — Ai, with its component perpendicular to the 
graphene's plane being B± = d x A 2 — d 2 Ai. 

Consider an auxiliary single-particle Hamiltonian H: 

H = mM + *7b7iA. (5) 

where M is a Hermitian 4x4 matrix. When m = 0, H = 
Ho, for each spin state. Also, if M 2 — 1 = {M, 7o7i} = 0, 

H 2 = D\ + B(<j z ®<j z )+m 2 . (6) 

This is the case if either 

At = Mi = a(I 2 <g> a z ) + b(a x ® a x ) + c(o y ® a x ) (7) 

with real a, b, c which satisfy a 2 + b 2 + c 2 = 1, or 

M = M 2 =a z ®a z . (8) 

In either case the eigenvalues of H 2 are at 2nB + rn 2 , 
with n = 0, 1, 2, .... For n > this immediately implies 
that eigenvalues of H itself are at ±v / 2n_B + to 2 , with 
the degeneracies of B/tt per unit area being the same as 
for m = 0. For n — 0, on the other hand, an elementary 
calculation gives that the eigenvalues of H, a) for any 
Mi, are at ±|to|, each with halved degeneracy of B/2ir 
(per unit area), and, b) for M 2 , are at to, still with the 
full degeneracy of B/tt. The invariance of the spectrum 
in the first case under rotations of the unit vector (a, b, c) 
is the consequence of the 'crural' SU (2) symmetry of Ho 
generated by {735, 73, 7s}, where 73 = a x <g> a y , 75 = a y <g> 
a y , and 735 = 17375, for example 0], Any specific 

choice of Mi in H breaks this SU{2) down to U(l), and 
leads to the same eigenvalues. M 2 preserves the chiral 
symmetry, and hence implies a different spectrum. 

Assuming that such an effective chiral-symmetry- 
breaking gap to becomes generated by interactions and 
then splitting the resulting energy levels further by the 
standard Zeeman effect could thus lead to the QH states 
at all even/, / = 0, and / = ±1, in accordance with ex- 
periment @]. Here the filling factor is / = 2ttN/B, with 
N as the number of electrons measured from a half-filled 
band. This still, however, leaves a complete freedom of 
choice of the vector (a,b,c). In particular, this choice 
may differ for the two spin states as well. For example, 
at B = and at strong V, one expects the interactions 
to prefer a = 1 for both projections of spin, i. e. the 
CDW 0. For strong U and at B = 0, on the other 
hand, it is a = a, i. e. the AF, that has lower energy 
[Tj. a = would correspond to the competing 'Kekule' 
ordering [l6j]. Which direction on the chiral manifold is 
actually chosen by the system is thus obviously a ques- 
tion of dynamics, to which we turn next. 

The low energy Lagrangian for the extended Hubbard 
model may be written as [7| 

L = $5^* <r 7 p D^-53(r + a-5,fl)*t* a (9) 
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where = ^\{x,t)^q : Dq = —id T , with r as the imag- 
inary time, fj, = 0, 1, 2, and 70 = I2 <£> 0^ , 71 = er z ® tjj,, 
and 72 = I2 ® (7 X . Here g a w C/a 2 /8, ,g c w (3V^ - U)a 2 /8, 
g z is the (dimensionless) effective g-factor of the electron, 
and B = {B']_ + B' 2 ) 1 / 2 , with B\\ as the field's in-plane 
component. We have retained only the two least irrele- 
vant short-range interactions among those present in the 
full effective Lagrangian at B = Q. This will prove 
justified at the magnetic fields of interest, as discussed 
shortly. We have also set the mass of the electron m e to 
unity. 

Performing the Hubbard-Stratonovich transformation 
and neglecting the quantum fluctuations the free energy 
per unit area and at T = may be written as 

2 2 

F(m c ,m a )-F(0,0) = ^ + + (10) 

where m a = m c + erm a , and the function 

f tT (s,m) = 9(\m\-\r + ag z B\) + (11) 
C(sA 2 )(coth(sB_L) - 1). 

m. c /g c = J^a^^o)' m a/9a = E^ "^^*^)' are the 

CDW and the AF order parameters respectively. 
C(x) is the cutoff function that satisfies C(x — > 00) = 1 
and C(x — > 0) = 0. Summing over the LLs below a sharp 
cutoff in energy, for example, yields C{x) = 1 — e~ x . The 
first term in Eq. (11) represents the crucial zeroth LL 
contribution, and the second the remaining LLs. Zeeman 
energy is taken to always be smaller than the separation 
between the zeroth and the first LL, i. e. g z B < y/2B±_. 

Let us first consider the system at r = 0. The free 
energy is minimized by the solution of 

m+ m_ 4B±m+ f°° ds 2 , , in , 

-^=%rv ,n - /-(-.«»-). as) 

.9+ .9- 7r d/ ^ 7o s v 

where g± = g~ x ± g" 1 . Assume both g c and <? a to be 
weak, Ag c a <C 1, and positive, and g c > g a , for example. 
There are then three types of solutions: 

1) If m_|_ > m_ > .9^i?, then m+ = m_, i. e. m a = 
and m c = 2g c B±/w. This is the CDW. It exists when 
g c > Ttg z B/2B±. For g a > g c , of course, one finds m c = 
and m a = 2g a B^/it, i. e. the AF (Fig. 2). This, in 
particular, includes the case of the pure Hubbard model 
with V = 0. The linear dependence on B± reflects the 
proportionality of either m to the degeneracy of the LLs. 

2) For g z B > m + > m_ one finds a paramagnet with 
rn c = m a = 0, unless the stronger coupling exceeds the 



B = critical value of 7r/8A, which lies outside the weak- 
coupling regime. 

3) Finally, when m+ > g z B > the solution is 

Tn x — g x B±/ir, where x = c, a. Both the CDW and AF 
order parameters are finite, which would corresponds to a 
ferrimagnet. This solution exists for g c + g a > wgzB / B±, 
and \g c - g a \ < irg z B/B ± . 

At weak coupling only the first term, representing the 
zeroth LL in Eq. (11), actually matters for the gap m. 
The second term would become important only at strong 
couplings, or in the limit B± — ► 0. Comparing the ener- 
gies of the four possible solutions, we find only three sta- 
ble phases represented in Fig. 2. All three have the filling 
factor / = 0. Even if both g c and g a are positive it is 
always energetically favorable to open a single but larger 
gap that corresponds to the dominant coupling. The 
aforementioned ferrimagnet, which would have / = 1, is 
thus never the minimum of the energy at r — 0. This may 
also be understood as follows. At low magnetic fields, 
such that Ib S> a where Ib = l/y/B± is the magnetic 
length in our units, the flow of the couplings implied by 
the invariance of the gaps in Eqs. (12)-(13) with respect 
to change of the cutoff A, in the regime A> 1/Ib is the 
same as at B± =0 17]. All weak couplings are thus ir- 
relevant at intermediate length scales between the lattice 
constant and the magnetic length 0] . At B± 7^ the flow 
of the interaction couplings towards the stable Gaussian 
fixed point is terminated at the cutoff ~ 1/Ib , with the 
least irrelevant coupling being left as dominant in the 
low energy theory. This single surviving coupling then 
selects and generates the 'relativistic' gap at B± ^ 0. 
All other more irrelevant couplings can be neglected, as 
we partially did in Eq. (9), from the outset. Inciden- 
tally, this also justifies the Hartree approximation to the 
free energy utilized in Eq. (10): as long as the ground 
state is semimetallic and gapless at B± = 0, for a weak 
enough magnetic field the low energy theory at the cutoff 
~ 1/Ib *C 1/a is indeed weakly interacting, as assumed. 

It is useful to display explicitly the relevant energy 
scales in the problem. First, the 'relativistic' gap is 
to re (U/8)(B ± /B ), where B = 1/a 2 « 10 5 T is the 
characteristic scale for the magnetic field set by the lat- 
tice constant. Assuming U ~ lOe^ gives m « ImeV 
for B± ~ 10T. The LL separation, on the other hand, 
is Dl ~ t(Bj_/ Bo) 1 / 2 , and thus by roughly two orders 
of magnitude larger. The Zeeman energy, on the other 
hand, for g z ps 1 is D z ps (B/10 4 T)eV, and thus of a 
similar size as m. Since the widths of the plateaus at 
/ = 0, ±1,±4 are proportional to either m or D z , both 
much smaller than Dl, this would naturally explain why 
these QH states require higher magnetic fields to become 
discernible in presence of some fixed disorder-induced 
broadening of the LLs. 

The above discrepancy between the kinetic and the 
interaction energy scales is only the first in a general hi- 
erarchy of energy scales that is implied by the stability 
of the semimetallic fixed point at zero magnetic field. As 
discussed above, the magnetic field affects the flow of 
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the interaction couplings only at the length scales above 
the magnetic length I. The development of an order pa- 
rameter m x in some channel may be understood as the 
divergence of the corresponding coupling constant x at 
the length scale L x — l/m x ^> I. This is determined 
by x, via essentially dimensional analysis, by the relation 
m x l w xl~ n , where xl~ n is the dimensionless coupling, 
with l/l now playing the role of the ultraviolet cutoff. 
For x = g and n = 1 this way we reproduce the above 
relation for the interaction gap to, and for x = t and 
n = we get the characteristic kinetic energy scale Dj,. 
The crucial point is that more irrelevant interactions at 
-6 = 0, with higher negative dimensions n, translate at 
B / into energy scales that depend on higher pow- 
ers of B/B Q : m x ~ (B/BY n+1 ^ 2 . In particular, this is 
what justifies the replacement of a realistic finite-range 
interaction with the simpler S- function in Eq. (9). 

At r = there are thus two fundamentally different 
ground states: one at larger interactions that breaks the 
A — B sublattice symmetry, either in the charge (CDW) 
or the spin (AF) channel, with vanishing magnetization, 
and the other magnetic, at weaker interactions, with the 
full translational symmetry of the lattice. Both are in- 
compressible and yield the plateau in the Hall conduc- 
tivity at / = 0. If the parameters are such to place the 
system for Bu = into the former, increasing Bu would 
eventually cause the transition into the latter state. This 
certainly holds for the CDW, as well as for the AF with 
its Neel order directed along the magnetic field, as we 
assumed here jl 81 ] . 

Hereafter we will thus retain only the larger quartic 
coupling, call it g x , and its corresponding order param- 
eter in Eq. (10), to. At small m the free energy then 
becomes 

2 

AF=^(l + 0(Ag x )) + 0(m 4 ) (14) 

_^M(0(| TO | _ \ r + D A) + ( | TO | _ | r _ Dz \)) 

where D z = g z B. The non-analytic term ~ |m| comes 
from the zeroth LL. In the weak-coupling regime possible 
local minima are then at to/too — 0, 1/2, 1, with too = 
2g x B±/Tt. Determination of the global minimum of the 
free energy at a finite chemical potential then leads to 
the phase diagram presented in Fig. 1. Let us determine 
the filling factor when m = mo/2 solution is stable, i. e. 
in the strip r + D z > tuq and \r — D z \ < mo/2 in Fig. 1. 
From the first inequality we see that r + D z > m, and 
thus r > m — D z > — m — D z , so the values of two out 
of four spin- and sublattice-degeneracy-resolved energies 
of the zeroth LLs lie below the chemical potential. The 
second inequality implies that between the two remaining 
energies one is above and the other below the chemical 
potential. To see this, first assume r > D z . Then D z — 
m < r < D z + to, as announced. If r < D z , on the other 
hand, D z + m > r > D z — m again, so in either case 
/ = 1. The determination of the filling factor for the 
remaining states is even simpler and follows analogously. 



For / > 2 the 'relativistic' gap is m = 0, and there 
is a direct transition between / = to / = 2 at 
D z /mo < 1/4. For a fixed and larger D z , by increasing 
the chemical potential the system always passes through 
the intermediate / = 1 QH state, which is magnetic and 
breaks the discrete sublattice symmetry. The width of 
/ = 1 state is 2D Z — (mo/2) for D z /mo < 3/4, or too for 
Dj/fflo > 3/4. All transitions are discontinuous. As the 
ratio D z /mo can be changed by varying Bu, it follows 
that at a large enough Bu the width and the activation 
energy of the / = ±1 state becomes ~ B± and inde- 
pendent of -By , which should be experimentally testable. 
Since in the experiment Q the width of the plateau at 
/ = appears to be somewhat larger than at / = ±1 
and / = ±4, the latter being always 2D Z , we speculate 
that 1/4 < D z /mo < 1/2. If this is indeed the case the 
activation energy at / = should first decrease before 
increasing with Bu , whereas that of / = 1 would increase 
and saturate. 

The mechanism of 'magnetic catalysis' of the 'relativis- 
tic' gap Q utilized here has also been considered recently 
as an explanation of the integer QH effect in graphene in 
ref. [l^ . The crucial difference from the present work 
is that only ~ 1/r tail of the Coulomb interaction was 
included, which could produce only the CDW, and that 
with the large gap to ~ \/B± ~ D^. To have m much 
smaller from the LL separation, as observed experimen- 
tally, requires then what appears to be an unrealistically 
strong screening of the Coulomb interaction by the sub- 
strate [l2[ ■ In this theory the gap is also always inhibited 
by a finite chemical potential, and thus cannot exist at 
/ = 1 if it did not already at f = 0. This would imply 
that the activation gap in [l2| at / = 1 is always the 
Zeeman energy. Both results are in sharp contradiction 
with ours. 

The QH effect in graphene was also recently discussed 
in [l3j . The principal difference from the present work is 
that the AF ordering was entirely neglected, so that the 
existence of the QH state at / = 1 in [l3| is always due 
to the CDW, which then requires a large enough V. In 
contrast, in the present theory / = 1 state can exist even 
for V = 0, when it is due to the AF ordering. In particu- 
lar, setting V = g z = we find the pure (V = 0) weakly 
coupled Hubbard model at half filling and in magnetic 
field to be an AF, and not the Stoner's ferromagnet, as 
assumed in [13j. This is because a finite AF order pa- 
rameter, unlike the magnetization, besides splitting the 
zeroth LL, also lowers the energy of all other occupied 
LLs. Put differently, in considering the dynamics within 
only the zeroth LL, the couplings need to be renormal- 
ized up to the length scale of the magnetic length Is first. 
Then g a (lB) S> 9}{Ib) always, where 3/ is the coupling in 
the ferromagnetic channel pj], for this reason omitted in 
Eq. (9). Retaining such a ferromagnetic coupling may be 
shown to only increase the activation energy when to — 
to D z + 2g f Bj_/TT. 

CDW formation due to lattice distortion has been pro- 
posed as a mechanism behind the QH effect in graphene 
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in [14j , This explanation differs from the present and 
all the other ones in that Coulomb repulsion plays essen- 
tially no role in it. In particular, an / = 1 incompressible 
state would in that scenario appear with increase of the 
chemical potential for any Zeeman energy in contrast to 
our Fig. 1, where for weak Zeeman energy there is a 
direct transition from / = to / = 2. 

Finally, the present mechanism differs essentially from 
the recent proposal [15] , in which disorder is invoked to 
explain the absence of QH states at odd / ^ ±1. 

In this paper we have neglected entirely the effect of 
disorder. As usual, it will broaden the LLs, and thus pro- 
vide an intrinsic energy scale that needs to be exceeded 
in order for an incompressible state to be resolved. We 
believe that this is why the states at / = 0, / = ±1, 
and / = ±4 become visible only at higher fields than the 
main sequence at / = 4n + 2, since their energy gaps 
according to our scenario are inherently smaller. They 
may be therefore understood as the fine structure of the 
QH effect in graphene. 

To summarize, we assumed that the principal effect of 
the Coulomb interaction in graphene is to introduce the 
on-site and nearest-neighbor repulsion for electrons on a 
honeycomb lattice. Postulating further a semi-metallic 



ground state in zero field, we argued that at a finite field 
there exists a hierarchy of energy scales, determined by 
the degrees of irrelevancy of the corresponding couplings 
at B = 0. In this work we considered only the effects 
of the least irrelevant interaction, which we expect to 
give the leading correction to the Landau level struc- 
ture for the non-interacting electrons. The phase dia- 
gram of graphene at laboratory magnetic fields ~ 10T is 
proposed. The theory predicts the incompressible states 
at all even integer filling factors including / = 0, and at 
/ = ±1. The ground state of the system at / = either 
breaks the sublattice symmetry, in either charge or spin 
channels, or is magnetic, depending on the magnitude 
of the Zeeman energy. At / = ±1 the system is always 
in the translational-symmctry-breaking phase and with 
finite magnetization, whereas at |/| > 2 the sublattice 
symmetry is preserved. The phase diagram in the Zee- 
man energy - chemical potential plane is proposed, and 
several of its features that may be used to distinguish the 
present from other proposed scenarios are discussed. 

The author is grateful to J. Alicea, J.-N. Fuchs. V. 
Gusynin, and S. Sharapov for useful discussions and cor- 
respondence. This work was supported by the NSERC 
of Canada. 
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